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We present strictly efficient scliemes for scalable measurement-based quantum computing using continuous- 
variable systems: These schemes are based on suitable non-Gaussian resource states, ones that can be prepared 
using interactions of light with matter systems or even purely optically. Merely Gaussian measurements such as 
optical homodyning as well as photon counting measurements are required, on individual sites. These schemes 
overcome limitations posed by Gaussian cluster states, which are known not to be universal for quantum com- 
putations of unbounded length, unless one is willing to scale the degree of squeezing with the total system 
size. We establish a framework derived from tensor networks and matrix product states with infinite physical 
dimension and finite auxiliary dimension general enough to provide a framework for such schemes. Since in the 
discussed schemes the logical encoding is finite-dimensional, tools of error correction are applicable. We also 
identify some further limitations for any continuous-variable computing scheme from which one can argue that 
no substantially easier ways of continuous-variable measurement-based computing than the presented one can 
exist. 



INTRODUCTION 

To realize a quantum computer in the circuit model, it is 
crucial to have precise control over each of the carriers of 
quantum information. In addition to keeping the quantum 
state of the computer protected from unavoidable noise in- 
duced by the environment, it is necessary to implement suit- 
able quantum gates, usually in the form of one and two qubits 
gates. The latter is of particular difficulty, especially when 
optical quantum systems are used because their interaction is 
either weak or merely induced by measurements. 

The paradigm of measurement-based quantum computing 
(MBQC) as pioneered by Raussendorf and Briegel |[l] ID and 
substantially generalized by Gross and Eisert ||3] [H allows to 
get rid of the necessity of performing unitary operations to 
implement a quantum circuit. Instead, the actual computation 
is performed by preparing a multipartite entangled state, the 
resource, in a first step followed by adaptively chosen local 
measurements on this resource. The important improvement 
stems from the fact that the resource is universal, i.e., it can 
be prepared independently of the algorithm one wants to per- 
form. This means that the presumably difficult step, the one 
which involves entangling operations, can be performed off- 
line. This resource-preparation may also be probabilistic as it 
is possible to wait with the implementation of the algorithm 
until the resource is available. What is more, individual ad- 
dressing in the final read out step is anyway also required in 
the circuit model, so this step does not add any further diffi- 
culty to the scheme. 

Quantum computing based on continuous-variable (CV) 
quantum optical systems differs from the more conventional 
notion based on their discrete analogues in not making use of 
single photons as the caiTiers of quantum information ||5l|6l. 
Instead, it relies to a large extent on Gaussian states and their 
manipulation. The notable advantage is that Gaussian states 
are easier to prepare in the laboratory, the coiTesponding inter- 
actions are often stronger and easier to accomplish, and some 
measurements, e.g., homodyne detection, can be performed 



with an efficiency substantially surpassing the one of single- 
photon measurements. However, quantum information pro- 
tocols using Gaussian states only, Gaussian operations, and 
Gaussian measurements suffer from serious limitations as in 
this setting neither entanglement distillation |7-9| nor error 
correction against Gaussian errors ifTol is possible. Also, since 
one can easily efficiently keep track of first and second mo- 
ments, any Gaussian evolution of Gaussian states can be effi- 
ciently simulated on a classical computer 1 11 1, clearly ruling 
out as the possibility of universal quantum computing. 

MBQC based on Gaussian resource states has been exten- 
sively discussed in the literature [12-16| due to several fea- 
tures: They are a direct generalization of the well known qubit 
cluster-state to the continuous-variable regime, they allow for 
universal quantum computing with Gaussian measurements 
and a single non-Gaussian one, and they can be prepared 
with present day experimental techniques, as already demon- 
strated. However, if they are not formed out of idealized in- 
finitely squeezed states (ones that are not contained in Hilbert 
space and would require infinite energy in preparation) but 
rather from physical states possessing only finite squeezing, 
they suffer from exponentially decaying localizable entangle- 
ment. This limitation, which applies to any Gaussian resource 
state iiTespectively of the permitted class of measurements 
ifTTl . implies that full error correction and the machinery of 
fault-tolerant quantum computing fTSl [T9< . in parts yet to be 
developed for such continuous-variable systems, appears to be 
necessary even when both state preparation and measurements 
are perfect. Notably, those restrictions only apply when trying 
to perform a quantum computation of an unbounded length. 
For any finite length, there exists, for any required accuracy, 
a physical cluster state such that any quantum operation up to 
this length is possible with this accuracy. Because the widely 
assumed superior power of quantum computers compared to 
classical ones manifests itself most prominently in the scal- 
ing of the runtime with the problem size, the situation of a 
quantum computation with unbounded length is in the focus 
of attention in this work. 
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In addition to schemes based on Gaussian cluster states, 
two more classes of schemes relying on CV-quantum optics 
have been proposed. The first one is based on superpositions 
of coherent states which are called Schrodinger cat states or, 
when they have very low amplitude, kitten states |20|. They 
can be, in approximation, created by subtracting photons from 
squeezed states and allow for universal quantum comput- 
ing by passive operations such as beam-splitters and photon- 
counting detection only. However, they also suffer from quite 
severe problems: the probabilistic nature of the quantum op- 
erations stemming from the use of non-overlapping basis- 
elements seems to be the most fundamental one, giving rise to 
significant overheads. The second approach combines the ad- 
vantageous features of both discrete and continuous-variable 
optics ETl I22I . In these schemes, quantum information is 
carried by qubits and single qubit operations are performed 
directly on them. Two-qubit operations, in contrast, are per- 
formed by letting both qubits interact with a strong CV-mode, 
the qubus. 

With this article, we pursue two different, but complemen- 
tary goals: On the one hand, we aim at clarifying the bound- 
ary between settings where CV-MBQC is possible and such 
situations where this is not the case. For this reason, we 
will develop a general framework capable of describing quan- 
tum computation in the measurement-based paradigm, regard- 
less of the dimension of the carriers of quantum information. 
Within this picture, we can identify some serious limitations 
giving rise to challenges that have to be overcome. On the 
other hand, we introduce a first strictly efficient scheme rely- 
ing only on a simple controlled rotation, which can be realized 
by an atom-light interaction of the Jaynes-Cummings type or 
purely optically by the Kerr effect, for the creation of the re- 
source. On the measurement side, we require beamsplitters, 
phase-space displacers, and photon counting measurements. 
While this scheme is not fleshed out in all detail of its con- 
crete physical implementation, it should be clear that quantum 
optical implementations of such ideas are conceivable. 

The present article is organized as follows: First, we dis- 
cuss what operations are possible with continuous-variable 
quantum systems and introduce a framework based on ma- 
trix product states (MPS) to describe MBQC in a general set- 
ting. After a discussion of the properties which a CV-MBQC 
scheme needs to posses in order to be called theoretically 
efficient, we show the problems of achieving those require- 
ments with Gaussian measurements on non-Gaussian resource 
states. Last, we provide an example for a feasible scheme and 
discuss in detail how efficient MBQC can be performed in this 
situation. 



FEASIBLE PRIMITIVES FOR CV QUANTUM COMPUTING 

When performing tasks of quantum information with 
continuous-variable quantum optics different classes of oper- 
ations are considered to be of different difficulty. This is also 
true for measurements where the achievable efficiency greatly 



differs between the various methods. 



Gaussian operations 

Before continuing with the discussion, we remind the 
reader of some basic properties of Gaussian states and oper- 
ations while also taking the opportunity to set the notation. 
We consider a single light mode. The energy eigenvectors 
of the unit oscillator are denoted by |n) with n — 0, 1, . . .. 
The annihilation operator d acts on them according to d\n) — 
^/n|n — 1). The commutator relation with its adjoint, the cre- 
ation operator, is [a, d^ ~ 1, setting h = 1. The eigenvectors 
of the photon-number operator h = d^ci with n\k) = k\k) are 
called the number states or Fock states. One can now define 
the canonical operators or quadratures as 



(1) 



(2) 



and for 6 E [0, 2tt] the family of rotated quadrature operators 

qg = qcos{9) +psm{6). (3) 

The latter family of observables is the one that captures ho- 
modyne detection |23|. Gaussian unitaries are the ones which 
can be written SiS U — e'^^'''^^ where H contains no terms 
in higher than quadratic order in q and p (or, equivalently, in 
d and d^). There are three classes of Gaussian single-mode 
unitary operators into which all Gaussian unitary gates can be 
decomposed. The first ones are corresponding to the applica- 
tion of the displacement operator D{a) = exp(ad^ — a*d) for 
a e (D. Such a transformation is reflected in the Heisenberg 
picture by a map of the form 

q^ q + V2Rea, (4) 
pi-^p + V2lma. (5) 

Optically, such a transformation can be realized by mixing the 
mode with a second mode, which is in a strong coherent state, 
on a beam splitter in the limit of vanishing reflectivity. In the 
second one are the transformations generated by the clockwise 
rotation operator R{0) ~ exp{—i9h), which can be realized 
by a phase shifter, acting on the canonical operators as q i-> 
qe, p ^ qg+T^/2, while in the last one are those generated by 
the squeezing operator 



S{0 = exp 



-id' 



^t2 



(6) 



acting as q H- > e' g, p 1— )■ e~'^p. The single -mode Gaussian uni- 
tary operations form a (non-compact) group which we denote 
by Ug- To complete the set of operators, which are neces- 
sary to implement arbitrary multi-mode Gaussian operations, 
we introduce the (absorption-free) beam splitter acting on two 
modes 1 and 2 by 



B — exp 



t t 
-{d\d2 - did\) 



(7) 
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where t = cos(0/2) and r = sm{9/2) are the transmission 
and reflection coefficient, respectively. 

The Gaussian operations can be divided into two classes: 
Phase shifters and beam splitters do not change the total num- 
ber of photons and are called passive operations. To imple- 
ment a squeezer or a displacer, on the other hand, one does 
need additional photon sources which makes them more diffi- 
cult to realize. 

Other Hamiltonian building blocl(S 

There are also non-Gaussian operations which are within 
the realm of present experimental techniques: The cross Kerr 
effect gives rise to a non-linear coupling between two modes; 
its Hamiltonian reads 

H = Xni®n2. (8) 

It can be realized by transmitting two modes together through 
an optically non-linear medium. As non-linear optical effects 
are weak and absorption needs to be kept low, the achievable 
effective coupling strength is quite small, i.e., tx <C 1, where 
t is the interaction time. For the most relevant situation, where 
one of the modes carries a state with many photons while 
the other one is in the single photon regime, the Hamiltonian 
(j8]l has already been experimentally implemented to perform 
quantum non-demolition measurement (QND) ll24ll . The re- 
lated Hamiltonian 

77 = x|l)(l|®n (9) 

for a positive x describes the coupling of a two-level atom to 
a light mode in the dispersive limit of the Jaynes-Cummings 
model. In this situation, the effective coupling strength can be 
increased by placing the mode and the atom inside an optical 
cavity. A particularly strong interaction can be achieved with a 
technique known as electromagnetically induced transparency 
(BIT) which is routinely used to exchange quantum informa- 
tion between between light modes and atomic vapors and for 
which also experiments with single atoms exist li27ij28J . 

Measurements 

We now turn to the kind of measurements that we will con- 
sider feasible for the purposes of this work. The most impor- 
tant Gaussian measurement scheme is homodyne dectection 
which correponds to the observable qo in Eq. ([3]l. It is real- 
ized by combining the mode with a strong laser, called the 
local oscillator, in an interferometer, measuring the intensities 
on both output ports, and subtracting the results. Another im- 
portant type of Gaussian measurement is eight-port homodyn- 
ing corresponding to a direct measurement of the Q-function, 
i.e., Qp{a) = (al/oja) for a € (D, where 

|a)=e-H^/2 5^^|n) (10) 



are the state vectors of the non-orthogonal and overcomplete 
coherent states |l23l. 

Non-Gaussian measurements are in many instances more 
difficult to realize than Gaussian ones, and usually with sig- 
nificantly lower detection efficiencies. The most feasible 
instance of a non-Gaussian measurement reflects a single- 
photon detector which can only distinguish between the ab- 
sence of photons and the presence of one or more photons. 
The corresponding POVM-elements are |0) (0| and 1 - |0) (0|. 
A photon-number resolving detector is more challenging to 
implement but with time -multiplexing |25| or superconduct- 
ing nano- wires ll26l . it is possible to perform photon counting 
for the first few number states with reasonable efficiency. 



FRAMEWORK FOR CV-MBQC 

In this section, we introduce a general framework to de- 
scribe measurement-based quantum computing. A MBQC 
scheme consists two ingredients: First, a resource state and 
second a set of possible, or allowed local measurements. We 
start with a description of quantum wires which are used for 
single-qudit processing 1 29 1 and discuss their coupling to fully 
universal resources afterwards. 



Matrix product states 

The formalism of matrix product states, which was orig- 
inally introduced to describe a certain class of entangled 
one-dimensional many-body states lf30l [311 is extremely use- 
ful to capture the essentials of quantum computing in the 
measurement based model |4|. We start by describing one- 
dimensional continuous-variable quantum wires which can be 
used to carry a single qubit of quantum information. The tech- 
niques we develop are independent of the actual physical im- 
plementation but we will often describe them with terms of 
quantum optics and have such a system in mind. 

We consider a chain of L quantum systems with dimension 
dp which are called lattice sites. We also allow for the sit- 
uation dp = oo which describes CV light modes. We say 
that a state is physical if it has finite mean energy, i.e., for a 
single mode: Tr ifip) < oo. We choose a basis (or a count- 
ably infinite Hilbert-space basis) {\i) : i — 1, . . . , dp} of (C^p, 
which we call the computational basis, associate to any basis 
element a _D-dimensional matrix A[i] e A4d{'C), and write a 
translationally invariant matrix product state (MPS) as 

dp 

{L\A[tL]...A[l2Wi]\R)\zL,...,l2,^l) 

IL,.... 12,11=1 

(11) 

where \L) ,\R) E (C^ . This vector space is the one where the 
quantum computation will take place and is called correlation 
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space. The matrices must fulfill the completeness relation 



(12) 



for rendering deterministic computation feasible. 

An important fact is that all physical states can be approx- 
imated arbitrary well by a finite-dimensional MPS. As any 
finite-dimensional state can be written as a MPS 131], it only 
remains to show that one can truncate states with finite en- 
ergy. We write the mean total photon number of a state p with 
k modes as 



C30 

-/Vmoan = ^ ^^"1 
ni,...,nfc=0 



■ • + nk)pni,...,nk;ni,...,ni, 

oo 

^^(^max ~^ 1) ^ ^ Pni ,. ...nfc ;ni ,. ...nfc 5 



ni,...,nfc=n„ax + l 



(13) 



where the matrix elements of p are the ones of the Fock basis. 

Let ptrunk be the non-normalized state obtained from p by 
setting Pni,....nk;n[,....n'^ = if One of the indices is larger 
than rimax- This state fulfills 



Ptrunc J _ 



k{r 



1) 



(14) 



Denoting the trace norm, which is the relevant norm for the 
distinguishability of quantum states, by || • we get from 
Ref. 1,32,1 



llPtri: 



Pill < 3 



k{n„ 



1 



1/2 



(15) 



If we assume the mean photon number (the energy) per 
mode to be constant, the truncation error is independent of 
the system size. This means, for a given required accuracy, 
we can choose a cut-off number rimax independently of the 
length of the intended computation. 



MPS as quantum wires 

We now show how single-qudit MBQC can be performed 
with such a matrix product state, following and extending the 
line of reasoning used in Ref. (4]. Let us assume, the first 
mode is measured in the computational basis with result k. 
Then, the remaining system is described by the state vector 

dp 

|«')(x J2 {L\A[iL].-.A[z2]A[k]\R)\iL,...,i2). m 

iL,...,i2 = l 

This can be interpreted as the action of the matrix A[k] on 
the right boundary vector \B). To interpret this as the action 
of a quantum gate it is necessary that A[k] is proportional 
to a unitary matrix, i.e., ^[fc]^^yl[fc] = p{k)l where p{k) is 



the probability with which the measurement result k is ob- 
tained. When the measurement basis is continuous, we denote 
by p the probability density while P denote the corresponding 
probability measure. Because we will also discuss measure- 
ments where the corresponding matrices are not proportional 
to unitary ones, we note that the probability in this more gen- 
eral situation is p{k) = \\A[k]\il])\\'^ , where j-)/)) is the state 
vector of the correlation system. If the measurement is not 
performed in the computational basis but in another one and 
the result corresponds to a projection to the state vector \x), 
the matrix applied on \B) reads 



Ab[x] = ^{x\i)A[i 



(17) 



Note that the basis B = {\x)} may be continuous and/or over- 
complete and that measurements in such bases naturally occur 
in CV-quantum optics. We also allow this basis to consist of 
improper eigenstates reflecting an idealized homodyning de- 
tection. We can now already note an almost trivial but very 
important necessary requirement for a single-qudit MBQC 
scheme: If there exists no allowed measurement basis such 
that A[x] is proportional to a unitary matrix for almost all x, 
it is not even possible to transport a single Z?-dimensional qu- 
dit. An equivalent formulation is that p{x) must not depend 
on the state vector \tp) of the correlation space. If this was 
the case, a measurement would yield information about ji/)), 
which would clearly destroy coherence. 



Sequential preparation 

A possible way of preparing a matrix product state with 
bond dimension D is a sequential interaction with a D- 
dimensional auxiliary system 1(311 [33l . This picture proves 
to be extraordinarily useful when discussing quantum wires. 
Let the interaction between this auxiliary system and a local 
physical system be described by 



U : (C^ (g) (C'" ~¥ (C^ (E) (C^" 



(18) 



and assume the physical system to be initialized in the state 
vector 1 5*). Then, the matrix elements of the MPS matrices 
read 



A[z],,, = i{k\®{i\)Ui\j)(g>\^))^ 



(19) 



In this picture, which is sketched in Figure [T] the correlation 
system is identified with the auxiliary system, and a single 
qudit gate is performed by first letting the auxiliary system in- 
teract with a mode in a known state and then measuring the 
mode. It is important to note that this picture works regard- 
less of the actual way the resource state is prepared. For this 
reason, we will use the picture of sequential preparation to 
be more intuitive while maintaining full generality. For ex- 
ample, resource states can also arise as ground states of local 
Hamiltonians or be prepared by the action of translationally 
invariant quantum cellular automata 



Encoding of quantum information in correlation spaces 



If the goal is not to process a Z?-dimensional qudit but only 
a d-dimensional one with d < D, the above requirement is too 
strong. The c?-dimensional space is called the computational 
subspace. We introduce the encoding 

V ■.(C'^-^(C" (20) 

with V'^V = Id and the encoded matrices 

B[T]=V^A[i]V , BeMd{<C). (21) 

If they are proportional to unitary matrices for all i, it is possi- 
ble to process a single d dimensional qudit. One can observe 
that 

{cj,\V^At,[iL]-.-AtAii]V\^lj} 

^{<f>\Bb,[iL]...BbAHm (22) 

where \tp) € (C^. The subscripts bi denote the chosen 
measurement basis. This means, for all measurements, the 
state with the matrices B behaves exactly like the one with the 
matrices A. Thus, we can just proceed as if we had a state with 
a d-dimensional correlation space from the very beginning. 
This situation is not the most general one. The isometry 



(20 1 can depend on the measurement basis and the outcome 
of the previous steps. In this case, the computational subspace 
is not fixed but changes during the computation. It is reason- 
able to demand that at least the initial and the final encoding 
coincide so that it is possible to work with a fixed in- and out- 



coupling, provided by V. In this case Eq. (21 1 becomes 



mUt, [tL]Vt,_, [iL-i] . . . [ii]Vl [zi] A, [^l]V\iJ) 

= {cb\BtAiL] ■ ■ ■ BbAnM) (23) 

where B[ik] = V^^ [ik-i]- When several quan- 

tum wires are used to perform a quantum computation, it is 
also necessary to demand the changing encoding to return 
to the initial encoding given by V, whenever a coupling oc- 
curs. If this was not the case, the couplings would have to 
depend on the history of measurements and their results in 
both of the wires. This would clearly be against the spirit of a 
measurement-based protocol, where all the adaptation lies in 
the choice of measurement bases only. 

The infinite-dimensional matrix product states used in this 
work bear some resemblance to the continuous matrix prod- 
uct states which where introduced in Refs. ll36l l37l to de- 
scribe one-dimensional bosonic quantum fields. In contrast to 
the discrete structure of bosonic modes described here, they 
are formulated without an underlying lattice. However, such 
cMPS can be obtained by a suitable continuum limit from 
infinite-dimensional matrix product states. In this limit, only 
two independent MPS matrices remain while the rest can be 
calculated from them which is an important difference to the 
present situation. 




FIG. 1. Interpretation of single-qudit MBQC as sequential interac- 
tion with an ancilla. The physical sites are initially in the state vector 
1^) and interact with the correlation system through U. Depending 
on the measurement results in Mi and M2, the unitary gate U is 
applied on | . 



Coupling of wires 

To go beyond single-qudit processing and achieve quantum 
computational universality, wires have to be coupled. Physi- 
cally, this can be performed in different ways. One possibility 
is to let the auxiliary systems of the sequential preparation 
scheme interact. However, this might be very challenging to 
do in reality, both when working with atoms, due to the diffi- 
culty of controlling them, and with light, due to weak non- 
linear interaction. For this reason, joint measurements are 
preferable to couple two wires which is also shown in Fig- 
ure |2] In an optical situation, this could mean to combine 
two modes belonging to two different wires on a beam split- 
ter and measuring both output ports. This is a broadening of 
the usual definition of MBQC where only local measurements 
are performed while we now also allow for two-local ones. 
On the other hand, there also exists a closely related third way 
which is also based on the beam-splitter interaction but strictly 
stays within the measurement-based paradigm. The idea is to 
perform the coupling independently of the executed algorithm 
according to some fixed scheme. This step then belongs to the 
preparation of the resource while the subsequently performed 
measurements are purely local. 

Let W be the coupling between the two modes, and 
{12:2)} the measurement bases after interaction. Then, the ap- 
plied operation reads 



A[: 



Xi,X2 = 



[xi,x-2\W\i,3)A\i\® A{f^ 



(24) 



REQUIREMENTS FOR EFFICIENT MBQC SCHEMES 

In this Section, we give three condition a MBQC-scheme 
has to fulfill to be called efficient and comment on their im- 
portance. To facilitate the derivation, we restrict ourselves to 
the, arguably most important, case of D = 2 while stressing 
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FIG. 2. Coupling of two quantum wires. Tlie interaction of two 
physical systems through W and the subsequent measurements (Mi 
and M2) induces a entangling gate on the two correlation systems. 



be efficiently approximated. Because it depends on the mea- 
surement outcome which operation is applied, we cannot hope 
for this to be possible deterministically. Instead, the measure- 
ments induce a random walk on S1LJ2 which can be controlled 
by choosing the measurement basis in the next step depending 
on the gate implemented so far. For efficient quantum com- 
puting, the average number of steps must not increase too fast 
with the desired accuracy. We set 



where B denotes the measurement basis and define 

n 



(26) 



(27) 



that the requirements for universal computing on qudits with 
larger dimension are completely analogous. When discussing 
requirements, one has to distinguish between two different no- 
tions: One the one hand, there are theoretical requirements, 
of the kind that, if they are met, efficient scalable quantum 
computing is possible in the absence of noise. Going away 
from this idealized situation, on the other hand, two new fea- 
tures arise. Firstly, in practice, any quantum operation is af- 
fected by noise, which error-correction may allow to counter 
to some extent. Secondly, any actual quantum computation is, 
needless to say, finite. Thus, the practical requirements that 
might be sufficient to sustain a computation may differ from 
the theoretical ones. However, as we are primarily interested 
in scalable quantum computing and for reasons of clarity we 
will focus on the theoretical requirements and comment on 
other practical implications whenever it is in order. 



Conditions 

We start with the theoretical requirements for quantum 
wires (which process single qubits) and state the first neces- 
sary condition. 

Requirement 1 (Transport) Let !-(/;„) be the state vector of 
the correlation system after n steps in the quantum wire when 
the initial state vector was \'ip). We say that transport of a 
single qubit is possible if there exists some e > such that for 
all n 



F ■ — 



nun •- max Tti\\\\{lp\U\%l)„)Y 

U&V2 W 



> 



1 



(25) 



Here, U is the best attempt at undoing the action of the quan- 
tum wire on This requirement ensures that computations 
are not limited in length. If encodings of single logical qubits 
in multiple wires are used, this requirement can be weakened 
to 1/e — poly(n). 

For universal single-qubit quantum computing to be possi- 
ble, it is not enough to merely transport quantum information. 
It is also necessary that any one-qubit unitary U G SU2 can 



Requirement 2 (Single qubit universality) Let U e §U2 

and e > 0. There exits a sequence (possibly adaptive) of 
bases Bi such that the expectation value for an approximation 
U = Ajs [k] fulfilling 



\\U -U\\ < e 
after at most n measurements in the wire satisfies 



E (71) = O polylog 



(28) 



(29) 



Here, 1 1 • 1 1 denotes the operator norm, which is the meaningful 
figure of merit in this case. Of course, single qubit operations 
are not sufficient for universaUty but an entangling operation 
is also needed: 



Requirement 3 (Coupling) Denote with 
Cz = diag(l, 1,1,-1) 



(30) 



the controlled- Z gate. The expectation value for an approxi- 
mation C z fulfilling 



\Cz-C,\\<e 



(31) 



after at most n measurements on both of the quantum wires 
involved, where both single-site and two-site measurements 
are allowed, must satisfy 



E (n) = O polylog 



(32) 



Requirements |2] and [3] guarantee that every quantum circuit 
can be simulated in correlation space with a polylogaiithmic 
overhead. This is the commonly used requirement in quantum 
computing. When one is content with a polynomial overhead, 
(|29l) and (|32l) can both be relaxed to 



= O poly 



(33) 
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for some constant c > 0. A computational model fulfilling 
only this weaker condition is still of interest when the ultimate 
aim is to perform a quantum algorithm which provides ex- 
ponential speed-up over the corresponding classical one (like 
Shor's algorithm). However, whenever the speed-up is only 
polynomial (like in Grover's search algorithm), the stronger 
polylog scaling is clearly necessary for the quantum algorithm 
to have any advantage. When one aims at the implementation 
of constant size quantum circuits only. Requirements 2 and 3 
are not necessary as one does not need to implement the quan- 
tum gates to arbitrary accuracy but only to some e > which 
is determined by the total length of the computation. 

The last feature to demand from a MBQC model is the abil- 
ity to initialize the correlation system into some known state 
and to perform a measurement of the correlation space in the 
computational basis. 

Requirement 4 (Initialization and read-out) For efficient 
initialization of the correlation system to be possible, 
there must exist some measurement sequence such that, 
independently from its initial state vector 



IIIV'n)-|0)||<£. 



(34) 



Here \il)n) is the state vector after n steps in the wire and it 
must be true that E(n) = O (polylog 

For readout, it is necessary that after n steps with 



E (n) = O polylog 



the combined normalized action 

1 



(35) 



(36) 



]=n 



fulfills either \\AB[k]~\0){0\\\ < e or ||AB[k] - |1>(1||| < e. 

The first requirement just means that one can approach one 
of the computational basis states fast enough, while the sec- 
ond one ensures an efficient implementation of an approx- 
imative measurement in the computational basis. While 
measurement-based computational schemes that do not re- 
spect these requirements may be conceivable in principle, 
within the framework presented here they are very natural in- 
deed and necessary for universal quantum computing. 



Consequences 

In the previous section we have discussed the requirements 
a MBQC scheme has to fulfill in order to be called effi- 
cient. We now show some consequences of these require- 
ments which will help us both to find classes of schemes 
which cannot be efficient and such that are. This seems an 
important enterprise in order to arrive at "no-go results", so 
to identify serious limitations that have to be circumvented in 



one way or the other. For many important situations, Fmin in 
Requirement [T] vanishes exponentially in n as the following 
observation shows, which by virtue of the above statement 
implies that efficient quantum computing is not possible. 

Observation 5 (Impossibility of transport) Given a basis 
B, let C C B, with P{C) > 0. If for all x € C, A[x] U 
with U G SU2, transport is impossible when measuring in 
this basis. 



Proof: We calculate the fidelity 



fmin{x) ■■= max min 

u[x]esv2 W 



p{x) 



(37) 



Because P{C) > and /min (a;) < 1 for all x G C, there 
exists a non-empty compact interval I C C such that there 
is some S < 1 with /inin(a^) < S for all x G I. Using an 
argument from Ref. |17|, one can show that this implies the 
exponential decay of Fnun, showing that Requirement[T]is not 
satisfied. Thus, for transport to be possible, there must exist a 
basis such that almost all measurement outcomes correspond 
to some A[x] with j4[x]^A[a::] cx 1. For transport over a fixed 
length. Observation |5] does not apply because even if the fi- 
delity decays exponentially fast, it might be still large enough 
for the task under question. 

We show a simple sufficient condition, which will be used 
later when discussing examples for which Requirement |2] is 
true: 

Observation 6 (Condition for an efficient random walk) 

Assume that for any target gate U G EV2 there exists 
some basis B and some C (Z B for any P{C) > p{e) with 
l/p(e) = O (polylog (1/e)) where e is the desired accuracy. 
When 



\\AB[k]~U\\<s 
for all k E C, Requirement^is fulfilled. 

Proof: The expected number of tries can be bounded by 

00 

=0 I polylog 



(38) 



(39) 



The most important important situation covered by Observa- 
tion|6]is the one where the measurements have discrete results 
and for every U G SU2 a basis exists which contains at least a 
single element which fulfills the requirements of Observation 
|6] Sufficient conditions for initialization and read-out of the 
correlation system to be possible are provided by the follow- 
ing observation: 

Observation 7 (Sufficient conditions for Requirement |4| 

Let there exist some S > I and a basis B containing some 
subset C d B of "non-unitary measurements". Assume that 



A[x] cx U[x] 



si{x) 

S2(x) 



(40) 
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where 6 > is a parameter characterizing the interaction 
strength. Taking a coherent state vector |a) with a G R, 
a > as input, go = \/2acos{9), po — \/2asm{6), a q- 
quadrature measurement with result x applies in correlation 
space 



A[x] — diag (e 



ipa{x- 



for all X € C, where U[x\, V € SU2 and Si{x), S2{x) are the 
(not necessarily ordered) singular values of A[x\ which fulfill 
either Si/s2 > 6 or S2/S1 > In this case, efficient read-out 
and initialization of the correlation system is possible. 

Proof: The situation is most transparent when S2{x) = 0. 
In this case, a measurement result x corresponds to a projec- 
tion of the correlation space state to the state vector V\0). In 
the more general case the measurement is not projective but 
merely weak. By performing the unitary operation t/L/^^[a;], 
which is efficiently possible due to Requirement |2] and re- 
peating the measurement, a projective measurement in the ba- 
sis {T/|0), can be approximated. This approximation is 
efficient, due to the existence of a finite gap between si and S2 
occurring with finite probability. The independence of 1/ of a; 
is crucial: If this is not the case, the basis in which the cor- 
relation space measurement will occur cannot be fixed. This 
will lead to a destruction of quantum information resulting in 
a failure of the MBQC scheme. Initialization of the corre- 
lation system in the state vector |0) can be performed in the 
very same manner One just performs the measurement pro- 
cedure outlined above after which the correlation system is in 
[/[a;] |0) or [/[a;] 1 1), depending on the outcome. Applying now 

the gates U[x]'^ or XU[x]\ respectively, achieves the required with probabiUty p^{x) = exp(-(a; - qq 



-ipo{x- 



(43) 



with probability p{x) = exp(— (x — 90)^)/%/^- Thus, Re- 
quirement [T] is fulfilled. However, efficient single-qubit gates 
are impossible due to the lack of control. The only possible 
measurement (up to squeezing) for which A[x\ is unitary are 
phase-space displacements followed by q-quadrature homo- 
dyne detection as for all others, the projectors do not fulfill 



mx)\e^'a)\ = mx)\e'^'a)\ 



(44) 



which implies that the probabilities of the measurements de- 
pend on the state of the correlation system and, therefore, the 
applied matrix is not unitary. Writing the displacement as 

A = (A, + iAp)/V2, we get 



A 



diag(e'(fo+^-)(--^),e- 



»(po-Ap)(^-«i±^) 



(45) 

Ag)2)/V^. Up to 



initialization. 



LIMITATIONS 

When devising a MBQC scheme, there are several funda- 
mental limitations concerning the use of Gaussian states and 
Gaussian measurements. The first important fact is that all 
protocols involving only Gaussian measurements on Gaussian 
states can be simulated efficiently on a classical computer, rul- 
ing out the possibility for any quantum speed-up [11 1. How- 
ever, for an ideal Gaussian cluster state, universality can be 
achieved by adding a single non-Gaussian measurement to 
the toolbox |12]. On the other hand, any physical Gaussian 
quantum wire, including an one-dimensional Gaussian clus- 
ter state cannot fulfill Requirement [T] even when allowing for 
non-Gaussian measurements tl7l. 



a redefinition go '-^ 90 + A^ and some global phase, ( 45 1 is 
constant in A, i.e., there is no way of controlling the random 
walk, making it impossible to meet Requirement [2] The same 
argument holds for any controlled Gaussian operation and for 
any Gaussian input state. 

When allowing for non-Gaussian input states, some amount 
of control is indeed possible, as the following example shows: 



Consider the interaction of Eq. ( 42 1 together with a superpo- 
sition of two different photon numbers as an input, — 
(|0) + |2))/\/2. Choosing 6 = 7r/4, the post-interaction state 
reads 

a|0)-^(|0)-z|2)) + 6|l)-^(|0)+z|2)) (46) 

and the normalized matrices Aq and Ap for the q- and p- 
quadrature measurements are 



Aq[x] =diag(e-'*(^\e**(^)), 
Ap[x\ =diag(e"^(^\e-**(^)), 



Controlled Gaussian operations 

We now discuss further limitations and consider the situa- 
tion where the interaction ( fTSj ) can be written as 



|0)(0|®Go + |l)(l|®Gi 



(41) 



with 



and 



(x) — arctan 



2a;2-l 



V2 



up to local unitaries on the correlation system, where 
GojGi S Ug- This class contains states which allow for 
transport of a qubit in a wire by Gaussian measurements only. 
An example is given by 



p{x) = (4x4 - Ax^ + 3) exp(-a;2)/(4V^). 



(47) 
(48) 

(49) 

(50) 



U 



1 exp(-i6'n) + |1)(1| (g) exp(i6'n) (42) 



Thus, there exist two Gaussian measurements which allow to 
control the random walk but even this control is not enough to 
meet Requirement [2] Because this is a general feature of any 
MBQC-scheme relying solely on Gaussian measurements we 
will discuss it in detail. 
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General resource states with Gaussian measurements 

We consider an MPS with D = 2 with the continuous fam- 
ily of matrices Ag [x] with a; G R. As the projectors describing 
eight-port homodyning are not orthogonal, they will turn the 
state of the remaining system into a mixed one when measur- 
ing a mode, and we do not need to consider this measurement. 
Because all relevant Gaussian measurements can be expressed 
as the application of squeezing, rotating, and displacing, fol- 
lowed by a g-quadrature measurement, we can write 

^0.a,,Ap,aN = cos{e)Ag[Xx + Ag] + sin(6')Ap[a;/A + Ap] 

(51) 

where 



1 



/2tt 



dx A„ 



(52) 



and where we have omitted global phase factors. If single- 
qubit transport is possible in this wire, we can, without loss 
of generality, assume that Ag[x] := ^0,0.0,0 [2^] is proportional 
to a unitary matrix for all relevant x. As both squeezing by 
a finite A and shifting in phase space does merely result in a 
redefinition of x in Eq. (51 1, we can restrict ourselves to the 
discussion of 



Aelx] =cos(0)ylq[a;] + sm{9) Ap[x] 



J2 e'^"Vn(x)AH 



(53) 



n=0 



where ipnix) are the energy eigenfunctions of the harmonic 
oscillator. We assume that the physical dimension is finite 
and discuss the required changes due to an infinite physical 
dimension afterwards. 

We now argue that Requirement [2] cannot be fulfilled in the 
present situation: 

Observation 8 (Impossibility of control) Let Ag [x] be as in 
( |53| ) with finite dp. Then, there exist some constants C > 0, 
A > 0, and Eq > such that for all U € SU2, € [0, tt], and 
all < e < So 



lPei\\Ae[x]-U\\ < e) < Ce' 



(54) 



where denotes the probability for the random variable x 
which corresponds to a quadrature measurement with angle 



If this is true, the expected number of tries to implement any 
unitary is bounded as lE(ne) > C~^{l/e)^ for Q < e < £q, 
and Requirement [2] is not fulfilled. The proof is given in the 
appendix. 

This result, which can be summarized as "no quantum 
computation with Gaussian measurements only" is comple- 
mentary to the one reported in Ref. ifTTll which can, in turn, 
be summarized as "no quantum computing with Gaussian 
states only". These results apply only to the model where 



all operations are noiseless and error-correction is not per- 
formed. Together, this means, both non-Gaussian states and 
non-Gaussian measurements are simultaneously necessary for 
continuous variable MBQC. However, the two limitations are 
of very different nature. The reason for the former is that there 
exists not enough localizable entanglement which means that 
the quantum information gets destroyed along the wire and 
transport or teleportation is not possible. The root of the latter 
limitation is the insufficient possibility to control what quan- 
tum gate is performed. 



Infinite-dimensional physical systems 

Up to now, we have only considered resource-states where 
the physical dimension is finite. We demonstrate the differ- 
ence between states with and without support on a finite num- 
ber of Fock states with the help of an example where the re- 
source state is described by the matrices 



A[x] = ^{f{x)l+zg{x)Z) 



(55) 



with J\f = / da; + |(7(x)p) to ensure the normaliza- 

tion condition 



J dxA[x]^A[x] ^ 1. 



(56) 



Both / and g are chosen to be real, even functions. Thus, 
also their Fourier transforms /, g are real and the matrices 
corresponding to qg are 



Ae[x] 



1 



7\f 

fix) sin(0)) 



((/(a;)cos(0) 

i{gix) co8{e) +gix)sm{9))Z), (57) 



which are all unitary. If the resource state is finite- 
dimensional, /, g, /, and g are of the form 



h{x) — poly(a;) exp(— 



(58) 



and the probability for Ag [x] to differ too much from the de- 
sired unitary is too large. If we allow the resource state to 
be infinite dimensional, the situation changes. In particular, 
it is possible that f{x) = const for \x\ < c for some con- 
stant c > 0. This allows to find a resource state where a q- 
quadrature measurement implements some gate with a finite 
probability. This is a striking example for a qualitative differ- 
ence between a proper continuous-variable state and all of its 
finite dimensional truncations. 

One might now be tempted to use this new insight to de- 
vise a MBQC scheme which allows for theoretically efficient 
single-qubit processing with Gaussian measurements only, but 
this is impossible. This can be immediately seen by truncat- 



ing the state according to (15 1. Then, the above results imply 
that single qubit MBQC is impossible when the energy per 
mode stays constant when scaling the system. If, on the other 
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hand, the energy per mode is increased with the length of the 
computation, this no-go result does not hold any longer. 

The reasons for the existence of such a severe limitation lies 
in the strong notion of efficiency. If Eq. ( [32| in Requirement 
[2]is weakened to (33 1, this restriction is not present anymore, 
and it is possible to devise a MBQC-scheme based on Gaus- 
sian measurements on a non-Gaussian resource. To do so. 



one can take the resource defined by the matrices (43 i. After 
changing the variable and neglecting global phases, one gets 

A[x\ - S{-poqo)S{~2pox) (59) 

with p{x) ~ exp(— a;^)/ -^Tr where 

S{<i)) = diag(exp(-i(^/2),exp(i0/2)) 

exp(— i(/)/2)diag(l, exp(z0)) (60) 

is the phase gate. Assume that the target gate is S{(t>). Chose 
go and po such that = —qopo and, therefore, ^[0] = S{9). 
Using ||S'(a) — >S'(/3)|| < \a — f3\, one can bound the proba- 
bility that the distance of the actually implemented gate to the 
desired one is larger than e as 

¥{\\A[x] ~ S{0)\\ >e)< ¥{\2pox\ > e) < Cpos (61) 

where C is some constant. This means, it is possible to ap- 
proximate an arbitrary phase gate in 0(1 /e) steps. Changing 
the interaction unitary to contain an additional Hadamard gate 
on the correlation space and using the decomposition detailed 
below together with the composition law for errors, it is easy 
to show that any U E SU2 can be approximated in 0((l/e)^) 
steps. 

FEASIBLE, EFFICIENT CV-MBQC SCHEMES 

In this Section, we present a scheme which allows for effi- 
cient CV-MBQC. The limitations discussed above force us to 
use both non-Gaussian resource states and non-Gaussian mea- 
surements. Even though the experimental requirements for its 
realization are indeed challenging, the scheme does only use 
primitives which all have akeady been demonstrated to be fea- 
sible in proof-of-principle experiments. 

Single qubit operations 



18 I to be 



We choose the interaction unitary in Eq 

U = {H®1) (|0)(0| (»exp(-i6'n) + |1) (1| exp(j6'n)) 

(62) 

where 9 > is a parameter. We initialize the modes in a 
coherent state vector |a) with a > 0. Eq. (62i describes a 



controlled rotation in phase space. The class of measurements 
corresponding to unitary evolution is given by the displaced 
photon-counting measurements, i.e., projections onto the state 



vector \ x,n) := D{x)\n) where a- G R. With the help of ( 1 1, 
the applied gates read 



AJn] 



H 



-in(p(x) 









where 



arctan 



a sm{9) 
^acos{9) + X 

The corresponding probabilities are 



Px{n) 



(63) 



(64) 



(65) 



and fulfill, for every x, the normalization condition 

Any single-qubit unitary operation can be efficiently ap- 
proximated with those operations: Every U E S1LJ2 can be 
written as 

U = S{cbi)HS{cl>2)HS{<l>3) (66) 
H is the Hadamard gate, and (pi £ [0, 2tt]. Rewriting this as 

U - HSiO)HSi(j)i)HSi^2)HSih), (67) 



we have decomposed every gate into four applications of ( 63 1 

a sin(0) 



Inverting ( 64 1 yields 



tan((/i(a;)) 



acos(6'). 



(68) 



To implement S{4>), one only has to choose x such that 
—2ip{x)n = (j) for some small n where n — and n = 1 
suffice and the former is only needed for (p ~ 0. For any 
fixed a 7^ and 9^0, the maximally necessary shift \x\ is 
bounded and with (|65|) there exists a lower bound < po < 



Px (n) for all relevant x and n. Using ( 67 1, one can see that 
the probability to obtain any chosen U E E1V2 in four steps is 
lower bounded by p"^. Thus, the conditions for Observation [6] 
are fulfilled and efficient single-qubit processing is possible. 



Coupling of wires 

We now turn to the coupling of two wires. Let the two 
correlation space qubits be prepared in \tp)j — aj\0) + bj\l) 
for j = 1,2. Then, the state vectors after the interaction with 
the light mode reads 



b,\l)\e^'a) 



(69) 



Displacing now both modes by A^. — a cos (6*), rotating the 
second mode by an angle of 7r/2, and setting 7 = as'm{9), 
the joint state vector becomes 

aia2|0,0>|7)N7)+ai^2|0,l)|7>|-«7> 

+ 61 a2 1 1 , 0) I - 7) I Z7) + 5i 62 1 1 , 1) I - 7) I - «7) • (70) 
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A balanced beam-splitter transforms two coherent states as 



\am ^ 



V2 V2 



(71) 



Applying this to ( 70 1 yields 



aia2|0,0)|^(l + z))|^(l-z)) 
+ai&2|0,l)|^(l-z))|^(l + *)) 
+M2|1,0)|^(-1 + *))|^(-1-^)) 
+6i62|l,l)|^(-l-*)>|^(-l + *)). 



(72) 



The matrices corresponding to a photon-counting measure- 
ment in both modes with results ni and n2 respectively read 



A[ni, ni] =diag (e ■* ,e * ,e * ,e 

X diag (e , e * , e * , e 



(73) 



which is, up to a local unitary, a controlled phase gate, i.e.. 



^[ni,n2] = diag(e 



T(ni + 7n2) i7r(3ni+5,i2) . 

,e 



) diag (1, e 



X diag (1,1,1, e~^''(-"^+"^)). (74) 

If ni + is even, A[ni, 71,2] is not entangling, and the local 
operations can be efficiently undone as described above. If 
rii + n2 is odd, a C^-gate is implemented up to local correc- 
tion. As this situation occurs with finite probability for any 
choice of a and 0, this means that Requirement |3]is fulfilled. 

It remains to show that initialization and read-out are possi- 
ble. To measure in the computational basis, we shift the mode 
after the interaction by A = a{cos{9) — i sin(6')) which turns 
the joint state of qubit and mode to 



a|0)|0) +6|l)|2msin(6')). 



(75) 



Counting the photons leads to a matrix fulfilling the condi- 
tions of Requirement |4] Note that there is some asymme- 
try: If the counted number of photons is larger than one, the 
qubit is in 1 1). If the result is zero, one has to undo the gate 
i/diag (1, —i) and repeat the procedure. Due to the gap be- 
tween 1 and I (0|2iQ: sin(6')) |, this is efficiently possible. This 
concludes the proof that the proposed scheme is efficient. 

As already mentioned, the coupling scheme discussed so 
far is not completely in the MBQC-paradigm because it in- 
volves the adaptive coupling of sites. However, this is not 
necessary: Consider two wires which are coupled every fc + 1 
sites by the beam-splitter described by Eq. fTT] ). Because 
C| = 1, two consecutive couplings can be undone by appro- 
priately choosing local operations between them. From Eq. 



FIG. 3. Non-adaptive coupling of four quantum wires where the 
wiggly lines denote the probabilistic measurement-based Cz gate. 
By scaling the distance between the single couplings polylogarithmi- 
cally in the length of the computation, the total probability of failure 
can be made arbitrarily small. 



probability po > that the desired gate is realized and prob- 
ability of failure in at least one of the wires is upper bounded 
by 2(1 — pq)^/^. Let n be total number of gates, the total 
probability of failure is bounded by 



Pfailurc 



< 2n(l -po)'"^* < Cine 



— C2fc 



(76) 



(65 I it follows that for every block of four sites, there exists a 



for suitable constants ci and C2. Thus, for a fixed probabil- 
ity of success, k has only to grow as fc = 0(polylog(n)), 
which means that the need for two-site measurements can be 
removed for the price of a polylog overhead. For a quantum 
circuit consisting of more than two wires, as depicted in Fig- 
ure |3] every wire should be alternately coupled to its left and 
right neighbor When the circuits consists of m quantum wire 
with all have a length n, the respective overhead behaves as 
fc — 0(polylog(nr7i)). 



Errors 

Even if a model fulfills the theoretical requirements it will 
still be affected by noise. Thus, error correction will be nec- 
essary to build a scalable quantum computer. In the scheme 
presented above, one plausible dominant source of errors is re- 
flecting the finite efficiency of the photon-counting measure- 
ment needed to perform both single and two qubit gates. If 
photon loss happens, a "wrong" operation is applied in the 
two-dimensional correlation space. When the rate of these 
qubit errors is low enough, techniques of error correction and 
fault-tolerant quantum computing may be applied to make a 
computation with unbounded length possible at the expense 
of a polylogarithmic overhead 1181 1191 . This contrasts with 
the errors stemming from the use of finitely squeezed states 
in the scheme based on Gaussian cluster states for which no 
method of error correction exists yet to achieve fault tolerance. 

Beside of the obvious errors stemming from decoherence 
and imperfections in the measurement procedure, the pre- 
sented continuous variable schemes have an important ad- 
ditional source of errors. For example, when performing a 
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shifted photon counting, uncertainty in the parameter of the 
shift operator D{x) results in uncertainty of the applied opera- 
tor We now investigate how those errors scale with the length 
of the computation. We restrict ourselves to the single-qubit 
case while noting that the error analysis for multiple qubits is 
completely analogous. Consider the situation where the tar- 
get gate reads U — 11"= i while the gate actually applied is 
A — n"=i -^i- We assume that the error induced by a single 
step is bounded, i.e., for all i, 

\\A-U,\\<e. (77) 

If all Ai are also unitary, the total error grows linearly in n and 
can be bounded as ||A — J7|| < ne. If this is not the case, an 
tight upper bound to the error is 

C/|| < n£(l + £)", (78) 

i.e., it grows, in the worst case, exponentially in n. In the 
MBQC-scheme based on Gaussian cluster states as discussed, 
e.g., in Refs. |12-15|, the implemented gates are not unitary 
for physical resource states which is the reason for the expo- 
nential decay of the transport fidelity ifTTl . 

CONCLUSION 

In this work, we have investigated the possibilities provided 
by measurement-based quantum computing with continuous- 
variable resources. After introducing a framework allowing 
for the description of a huge class of possible resource states, 
we have clarified which conditions are necessary for a CV- 
MBQC scheme to be viable. Those conditions lead to new 
limitations. Especially, we have shown that, without the use of 
not yet developed continuous-variable error correction, scal- 
able quantum computing with Gaussian measurements alone 
is impossible, even if the resource state is non-Gaussian, com- 
plementing the prior result which prohibits MBQC with Gaus- 
sian resources, even if the measurements are non-Gaussian. 
As the second main result, we gave an explicit example of 
an efficient MBQC scheme where the resource could be cre- 
ated by a simple interaction between light modes in coherent 
states and some qubit degrees of freedom. Processing is then 
performed by shifted photon counting while entangling two- 
qubit gates rely on simple beam splitters. An analysis of errors 
highlighted the qualitative differences of the two major kinds 
of errors. 

Two major conclusions can be drawn from our findings: 
First, continuous-variable measurement based quantum com- 
puting is an extraordinarily difficult enterprise. Second, 
when non-Gaussian resource states are combined with non- 
Gaussian measurements, efficient schemes do exist, and in 
this work, we introduce such first strictly efficient schemes. 
While this may be disappointing when aiming at a scalable 
optical quantum computer, one does need to keep in mind, 
that our statements about efficiency only address the asymp- 
totic behavior. Thus, the limitation do not at all rule out the 



usefulness of much simpler protocols in situations where only 
a few modes are used. Such situations include quantum re- 
peaters and hybrid entanglement-distribution protocols for ap- 
plications particularly in quantum cryptography. Hence, for 
such purposes, both the purely Gaussian setting as well as 
those non-Gaussian schemes that asymptotically do not give 
rise to universal quantum computing can well be feasible as 
elementary tools for entanglement distillation schemes. It is 
the hope that the framework established here fosters further 
such investigations. 
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Proof of Observation |8] 

To find an upper bound to P(||Ae[2;] — U\\ < e), we first 
lower bound the probability that |a;| > for some suitably 
chosen xq. Then, for |x| < xq, we upper bound ||^0[a;] — 
U\\ by a polynomial in x and use that the probability for a 
polynomial to be larger than e cannot grow too fast for small 
e. 

We calculate 



where IIOIU 



\\Ae[x]^U\\>\\ 



(79) 



X^i I'^ij P denotes the Frobenius norm. 



We now use d53|l and the fact that 



(80) 



where Hn{x) denotes the nth Hermite polynomial. We also 
use < y/x < a; + 1 and that P0{x) is upper-bounded by 
a constant independent of 9. This allows us to find non- 
negative polynomials fug and g{j,e, where the coefficients 
are continuous functions of U and 9 with — C/|| > 

fufi{x) / gif,o{x). As the probability density pe{x) is of the 
fom\pg{x) = poly(x) exp(— x^) we get 



> Ci log 



< e 



(81) 



for some suitable Ci > 0, C2 > 0. Now, we assume that 
\x\ < Ci log {C2/S). In this case, gjj g can be upper bounded 
by C3 (Ci log {C2/s)Y for some C3 > and some even r. 
Note that the maximal degree of all polynomials only depends 
on dp. Assume that fu,o has M zeros xi, . . . , xm- Around 
any of these zeros, fjj^e can be lower bounded by C/i{x — XiY 
for some C4 > and some even s. Using this and ( 8 1 1, we 
chose some Eq such that 

Ve{\\Ag[x\-U\\<e) <e (82) 
+ Yy<^ (^4(x - x,Y < C3 (Ci log (^^^ ^ . (83) 

for all e < Eq. Because P0{x) is upper bounded, there exist a 
constant C5 such that 



Fiix-x,y <e)< Cse 



l/s 



(84) 



for all e < Eq. Inserting (84i into (83i, one obtains p4\ . As 
all constants depend on U and 6* in a continuous way, and U 
and 9 are both taken from compact domains, there exist a set 
of constants such that j54| is simultaneously true for all U and 
9. Thus, ObservationlSlholds. 



